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Abstract 

The low-energy S'-wave component of the decay D + — * K~tt + tt + is studied by means of a 
chiral SU(3) x SU(3) effective theory. As far as the primary vertex is concerned, we allow for 
the possibility of either direct production of three pseudoscalar mesons or a meson and a scalar 
resonance. Special attention is paid to final state interactions associated with elastic meson-meson 
scattering. The corresponding two-body amplitude is unitarized by ressumming s-channel diagrams 
and can be expressed in terms of the usal phase shifts 5. This procedure preserves the chiral 
properties of the amplitude at low-energies. Final state interactions also involve another phase w, 
which describes intermediate two- meson propagation and is theoretically unambiguous. This phase 
is absent in the if-matrix approximation. Partial contributions to the decay amplitude involve a 
real term, another one with phase 5 and several others with phases 5 + lo. Our main result is a 
simple and almost model independent chiral generalization of the usual Breit-Wigner expression, 
suited to be used in analyses of production data involving scalar resonances. 
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I. INTRODUCTION 



Decays of heavy mesons have been recognized recently as important sources of information 
about scalar resonances, since the E791 experiment has produced solid evidence for a broad 
scalar-isoscalar state in D + — > 7r + 7r~7r + [1], know as the a. In the case of the reaction we are 
concerned with, namely D + — > K~tt + tt + , the E791 group included in their fit another scalar 
state, known as the k (or if* (800)), and concluded that this resonance was the dominant 
source of S'-wave n + K~ pairs[2]. This represented a turn point in our understanding of this 
reaction, which was hitherto thought to be dominated by the non-resonant background. The 
existence of the k was then confirmed in a different approach by FOCUS [3] and by the new 
high statistics results from CLEO [4]. Empirical information about the parameters of these 
resonances, especially the position of their poles in the complex energy plane, can only be 
obtained after analyses of large amounts of data organized into Dalitz plots. 

Low-energy mesons correspond to gentle deviations from the QCD vacuum and tend to 
be highly collective states, as indeed happens with pions or kaons, which are large objects. 
In the case of low-energy scalar resonances, one is entitled to expect that the influence of 
the vacuum will be a truly overwhelming one. This is probably one of the reasons why they 
are so broad and prove to be so elusive. 

Dalitz plots of .D-meson decays into three pseudoscalars usually contain several S- and 
P-wave resonances [1, 5], which couple by means of final state interactions (FSIs). This 
makes the disentangling of individual resonance properties both very involved and strongly 
dependent on particular forms of trial functions adopted. Theoretical ansatze employed in 
data interpretation usually rely on Breit-Wigner expressions, which proved to work well in 
many instances. A popular ansatz for the trial function has the form 

A = a e*>Ao+ £ a s n e^A s n + £ a p n e** A p n + ■ ■ • , (1) 

S— wave P— wave 

where the first term represents a non-resonant background and the amplitudes A e n are Breit- 
Wigner expressions for each resonance present in the final state. The masses and widths of 
well established states are used as input, whereas those of low-lying resonances, as well as the 
free parameters a\ and <p\, are fitted to data. As a consequence, these adjusted parameters 
acquire the status of empirical quantities. 

In the case of tctc scattering, the fact is well established that the lowest pole of the 
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amplitude is located roughly at ^/s (0.47 — i 0.29) GeV, whereas the S'-wave phase shift 
reaches ir/2 around v /s^7 ss 0.92 GeV. These findings are clearly at odds with the traditional 
use of Breit-Wigner expressions. An explanation for these seemingly paradoxical results was 
provided by Colangelo, Gasser and Leutwyler[6], who have shown that chiral symmetry 
requires a compromise between the polynomial nature of the amplitude at very low-energies 
and the vanishing of its real part at s n/2 , which is responsible for a large shift in the pole 
position. This kind of feature is inherent to the isoscalar channel. The resonance, which is 
a non-leading chiral effect, must always coexist with an important polynomial in s. In the 
framework of chiral symmetry, the first two terms in eq.(l), namely [a e l<t>0 Ao + af e 1 ^ Af], 
are not suited for describing low-energy interactions. The use of this kind of trial function 
is problematic and may give rise to results which are not reliable. 

The extraction of information from experiments involving scalar resonances must be 
performed in the best theoretical framework possible, as the quality of results for masses and 
coupling constants depend on the ansatze employed. In this work we propose an alternative 
form for the low-energy sector of the trial function, which can be used as a tool in analyses of 
the decay D + — > K~7t + tt + . Our motivation for choosing this particular reaction is two-fold. 
The first one is that scalar resonances in the final state occur just in the 7r + K~ subsystem, the 
number of possible couplings is relatively small, and the problem is simplified. The second 
is the availability of recent data analyses on this process, which could allow the testing of 
our results. Nevertheless, the general lines adopted here can be extended to other systems 
in a straightforward manner. Our theoretical model is based on standard SU (3) x SU (3) 
effective chiral lagrangians incorporating scalar resonances[7, 8]. 

We are also concerned with the presence of phases in theoretical ansatze, which are 
sometimes employed with no visible physical meaning. In the realm of two-body systems, 
the dynamical origin of phases is well understood. In particular, it is known that two-body 
rescattering gives rise to the elastic phase shift 5. When resonances are present, this phase 
is shared with the so called production amplitude, as dictated by Watson's theorem [9]. We 
show, in the sequence, that another phase is also present in the D + decay considered here, 
associated with two-meson intermediate states. 

The theoretical description of a heavy-meson decay into three pseudoscalars is neces- 
sarily complex. Nevertheless, we have made an effort to produce final expressions which 
incorporate a compromise between reliability and simplicity, so that they could be employed 
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directly in data analyses. Our paper is organized as follows. In section II, we review the 
basic lagrangians, which are used in the description of both the primary weak vertex and 
final state interactions. The former is discussed in section III, whereas the latter are dis- 
cussed in sections IV, V and VI. These results are assembled in section VII, where our main 
results are presented and analyzed. In section VIII, we display individual predictions from 
the various components of the decay amplitude in Dalitz plots, so as to produce a feeling of 
their dynamical content. Finally a summary and comprehensive conclusions are presented 
in section IX. Details concerning kinematical variables, the form of the two-body propagator 
and background interactions are left to appendices. 

II. DYNAMICS 

The reaction D + — > K~7r + 7r + involves both weak and strong processes. The former 
are associated with the isospin conserving quark transition c — > sW + , whereas the latter 
occur in final state interactions involving both pseudoscalar mesons and their resonances. 
In order to keep approximations under control, we remain, as much as possible, within the 
single theoretical framework provided by a chiral effective field theory. This choice is also 
motivated by the fact that we are concerned mostly with the low-energy sector of the trial 
function. 

• strong interactions: 

Meson-meson interactions are described by the chiral SU(3) x SU(3) lagrangian at 0(q 2 ) 
given by Gasser and Leutwyler[7], whereas couplings of scalar resonances to pseudoscalar 
mesons are taken from the work of Ecker, Gasser, Pich and De Rafael[8]. Keeping only 
relevant terms, one has 

£ {2) = ^ ( v M c/t wu + x i U + xU i) 

+ C d ( S Ufj, ) + C m ( S X+ ) 

+ c d Si ( m m ) + c m Si ( x+ ), (2) 

where ( • • • ) indicates the trace, F, Cd, c m , dd and c m are constants, S and Si represent 
scalar resonances and U is the pseudoscalar field. Using the definition u = U 1 ^ 2 , one has 

V»U = dpU, u^ = iv)V^Uu\ 

X = 2Ba, x+ =^ xu^ + ux^u . (3) 
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The field a incorporates the quark masses as external scalar sources and, in the isospin 
limit, is written in terms of the usual Gell-Mann matrices as a = a® I + a 8 with 
cr = (2m + m s )/3, cr 8 = (m — m s )/\/3 and m„ = m d = rh. We neglect 77-77' mixing and 
use Ba = (2M| + M 2 )/6 and B a 8 = -(M^- M 2 )/V3. The meson field is written as 
U = e % ®' F , $ = Xi (pi , and the leading order strong lagrangian becomes 

'4 



£(2) 



1 i? 

/iis fkis 4>i d^j <p k dfj,(pi + 
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24F 2 



+ °8 [ dkl8 + -dij8 5kl + 2dij m dklndg m n 



X <Pj <Pk <Pi 



+ 



2q 
F 2 



S'o d^4>i 



4c n 
F 2 



B [a 5ij + cr 8 d 8ij ] S < 



2 Cd c m 
+ d U j S b dpfa d^j - — B [3 cr daij + cr 8 (2 5 bi 5 j8 + 4 d bis d sjS )) S b <pi <pj , (4) 

where fak and d^ are the usual SU(3) constants. We use the conventions of ref. [10] for 
meson fields and R° = (Se+iS^)/ y/2 for the scalar resonance. In the sequence, the R° state 
is called k for simplicity and, in isospin space, one has 



|vr + /Y- ) = |3/2, 1/2 ) + y^2/3 |l/2, 1/2 ) ; 

Itt ^ ) = ^2/3 |3/2, 1/2 ) - y/l/3 [1/2, 1/2 ) . 



(5) 



weak interactions: 



w + 



K 



(a) 



D + 



(b) 



D + 




FIG. 1: (Color online) Weak amplitudes involving pseudoscalar mesons (dashed lines) and the 
scalar resonance (continuous line). 



Effective weak vertices contain a propagating W and are based on the processes shown 
in fig.l, which involve the following combinations of reactions: 
type (a): (D + -> tt KW + ) a and (W + -> 7r) a , 
type (b): (D+ -> X W + ) v and (W+ -> tt tt),, , 
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type(c): (D+ -> R W + ) a and (W+ -> 7r) a , 

where the labels t> and a refer to either vector or axial currents in couplings. In the want of 
a comprehensive theory, these weak vertices can be derived by means of appropriate hadronic 
currents and ad hoc phenomenological coupling constants. The latter can be estimated in 
semi-leptonic processes, by replacing the top n + in the figure by (£ + v^). Although feasible, 
the piecemeal implementation of this program is cumbersome, owing to the large number of 
possible isospin couplings and phase conventions. 

Fortunately a rather economic alternative is available, based on the group 377(4). Before 
proceeding, we would like to make clear that we are quite aware that SU (4) is not a good 
symmetry and therefore we are not advocating its use here. On the other hand, the inclusion 
of D mesons into the pseudoscalar multiplet does give rise, by means of the external source 
technique, to an effective lagrangian in which all currents and isospin couplings are generated 
automatically. We obtain W + couplings by using a SU (4) x SU (4) mesonic lagrangian for 
the weak sector and, at the end, break any commitments with a symmetry, by allowing 
phenomenological coupling constants into the Feynman rules. 

In the standard model lagrangian, weak couplings are given by C w = q 7 M [v^ + 75 a M ] q, 
where the external sources are written in terms of S77(4) matrices as t> M = v% A a /2 and 
a M = A a /2. The charged currents are 

*V + a M = 0, <- v,-a^ = -^= [wfVsuw + W-V^ , (6) 

where g is the weak coupling constant, W is the gauge field and V S ij{a) is the SU (4) sector of 
the Cabibbo-Kobayashi-Maskawa mixing matrix. Its explicit form, in terms of the Cabibbo 
angle 9c, reads 1 

V S u(4) = ^ [cos 9 C (Ai + + sin 9 C (A 4 + iA 5 ) 

- sin 9 C (An + i X 12 ) + cos 9 C (X 13 + % A M )] . (7) 

The weak effective lagrangian is obtained by replacing V M in eqs.(2) and (3) with V M U — 
d^U — i (Vft+a^U + i U{v^ — a M ). The term that interests is then given by 



2 4cd 
F S a i 9^1 + — fajs fkls 4>j <Pk d^i + — d ab i S b d^c 



,(8) 



where the factors within square brackets are respectively vector and axial hadronic currents. 



1 We use the SU{4) conventions of ref.[ll]. 
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III. WEAK VERTICES 

In this section we display the tree level amplitudes given in fig.l and include final state 
interactions afterwards. The weak vertices in which the D + participates involve the trans- 
formation of a c-quark into a s-quark and are Cabibbo allowed. As final state interactions 
allow the transition tt°K° — > n + K~, one also needs to consider vertices involving neutral 
mesons. The production of a single pion in — > n + (q') is represented by 

T- = i{Fg/2) cos^*?;. (9) 

Decays of the type (a) are based on the processes D + (P) — > K~(k) n + (q) 

and D + (P) -> W+ K°(k) ir°(q), given respectively by the amplitues ^/2/3 T^ K and 

-073 T£* K , with 

T^ K = -i (g/2 V6 F) cos 6 C (P - fc)„ ■ (10) 

In the case of vector couplings, one needs the vertices D + (P) — > W 7 "^ K°(k) and W 7 "^ — > 
7r + (g') 7r°(g), which read 

T™ = -{g/2y/2) cos9 c (P + k),, (11) 
7T = (ff/2) cos M<Z" 9% ■ (12) 

Finally, for the vertex D + (P) — > re, involving the scalar resonance, one has 

T™ = i(c d gV2/F) cos 0c ^- (13) 

Using = i (gV/Myy) for the W propagator and the definition GV = \pl g 2 /8M^ = 
1.166 x lCT 5 GeV 2 [12] , one finds the following weak amplitudes 
type (a): 

jD +( P ) ^ K~(k) n + (q) n + (q') : y/2/3 W a , (14) 
D + {P) - 7r°(g) 7r + (g') : -^/T/S W a , (15) 

W a = [<J a ] (Gf/a/3) cos 2 6 C (P-k)-q' , (16) 

type (b): 

jD +( P ) ^ ^0( fc ) ^ . W6 = _ [6b] Gp cog 2 ^ (p + k) . (g[ _ ql) ? (17) 



7 



type (c): 

D + (P) -> R(q s ) n + (q') : W c = - [S c ] 4 G F c d cos 2 6 C P ■ q\ (18) 

where ad /joc factors [Si] were introduced so to freeing results from any constraints imposed 
by S77(4) symmetry. 

IV. FSI: KERNEL 

Final state interactions are essential to structures observed in Dalitz plots, since they 
promote couplings among various channels and, in particular, give rise to widths of reso- 
nances. The final state considered here contains three mesons and a complete treatment 
of the problem is not possible. One is forced to employ approximations and we adopt the 
quasi-two-body approach, in which one of the final mesons acts as a mere spectator. As 
emerging pions have isospin 2 and no resonance is known in this channel, their interactions 
can be safely neglected and strong interactions are restricted to the ttK subsystem. 

These assumptions lead to the model given by fig. 2, in which all strong processes are 
incorporated into the amplitudes T and IT, representing respectively elastic scattering and 
production. As the isospin of the 7iK system can be either 1/2 or 3/2 and only the former 
couples with the k, one needs to consider three amplitudes, namely T 1/2 , T 3/2 and U 1/2 . The 
basic building blocks of both Tj and 11/ are kernels JCj, which describe elastic S-wave tcK 
scattering at tree level. These kernels are obtained by projecting out S'-wave components 
from tree-level amplitudes Tj, given by the diagrams of fig. 3. As far as chiral symmetry is 
concerned, the leading term is given by the 0(q 2 ) contact term whereas diagrams involving 
resonances are 0(q 4 ) corrections. 

We consider the process n(Q) K(K) — > ir(q) K{k) and the tree amplitudes are explicitly 
written in terms of the Mandelstam variables as 

T 1/2 = ^r 2 [4 S + 3t-4(M 2 + M|)] 

"ii 4 (8-MZ-M*) +Cm (5M 2 +4M|)/6] 2 

+ f°-f, 8 /6-f M 74, (19) 

f a/ 2 = - 2^2 i s ~ ( M - + M k)] + f t - f t/Q + T:/2 , (20) 
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FIG. 2: (Color online) Diagrams contributing to the decay D + — > i^~7r + 7r + ; (a) corresponds to a 
direct process, (b-e) involve the 7T.K" scattering amplitude T, and (f) depends on the production 
amplitude II. 
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FIG. 3: (Color online) Tree-level ttK amplitude; dashed and full lines represent respectively pseu- 
doscalar mesons and scalar resonances. 



(21) 



where 



rpQ 

1 t 



— l — ^ { [c d {t-2Ml) + 2 c m Ml] [c d (t-2M 2 K ) + 2 c m M 2 K ] } 



Tf = - 



t-ml F A 



{[c d (t-2Ml)-c m (2M|-llM2)/6] 



(22) 



x [c d (t-2M 2 K ) + c m (10M 2 K -M 2 )/6)} 
1 4 



— [c d (u-M^-M' K )+c m (AM' K + 5M^)/QY 



u—m 2 . F A 

The projection into S- waves is performed using results from appendix A and one finds 
1 



(23) 
(24) 



!/2 4_p2 



[(4-3p 2 /2) S -4 {Ml + Ml)] 



— [c d {ml-Ml-M 2 K )+c m {hMl+AMl) /6] 2 + B 1/2 , 



s — m 2 . F 4 
1 



(25) 
(26) 



with p = a/1 — 2 (M|- + M^)/s + (M^ — M 2 ) 2 / s 2 . The functions 5/ are smooth backgrounds 
given explicitly in appendix C. As discussed there, all t and u channel are small and can be 
either treated as four-point contact interactions or neglected. This gives rise to the effective 
structures shown in fig. 4, where contact terms in the kernels include both leading chiral 
contributions and non-resonant backgrounds. 



FIG. 4: (Color online) Effective structures of the kernels fCf, the cross hatched bubbles represent 
effective contact interactions, which include both the leading 0(q 2 ) contribution and background 
terms. 

In the isospin 1/2 channel, it is convenient to emphasize the role of the resonance by 
factorizing the s-channel denominator and writing[13] 

,2 



/C 



7 



1/2 



s—mt 



! = (3/F 4 )[c d ( m l-Ml-M 2 K ) + c m (5 M 2 +4M|) /6] 2 
{(1/4F 2 ) [(4 - 3p 2 /2) a - 4 (M 2 + M 2 K )] + B y2 } ( S -m 2 ) 



(27) 



Of course, in spite of differences in form, eqs.(25) and (27) have exactly the same content. 

As discussed in the introduction, we are interested in mapping low-energy degrees of 
freedom of the amplitude D + — ► K + 7t + tt + . This means that masses and coupling constants 
must be kept free, so that their values can be extracted from experiment. On the other hand, 
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in discussing qualitative features of our results, we need to fix somehow these free parameters. 
In this case, we choose: m K = 1.2 GeV, (q, c m ) = (3.2, 4.2) x 10~ 2 GeV, (q, c m ) = (1.8, 2.4) x 
10~ 2 GeV [8]. It is worth stressing that we are by no means recommending these values. 
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FIG. 5: (Color online) Kernels fZj (continuous lines) and the leading 0(q 2 ) contact chiral contri- 
bution(dashed lines). 



In fig. 5 we display the full kernels /C/, together with their leading 0(q 2 ) chiral components, 
and it is possible to note an important isopin dependence of the results. The dynamical 
structure of /C 1/2 involves three different regimes. At low energies, for s between threshold 
and ~ 0.6 GeV 2 , it is determined by chiral constraints whereas, as s increases, it becomes 
dominated by the first resonance pole. For larger values of s, effects associated with other 
resonances, not considered in this work, do show up. Therefore, with the choice m K = 
1.2 GeV, the upper limit of validity for our results is assumed to be y/s ~ 1.3 GeV. The 
kernel /C 3/2 is repulsive and monotonic. 



V. FSI: SCATTERING AMPLITUDE 



The elastic ttK scattering amplitudes Tj are derived by using the two-body irreducible 
kernels JCj into the Bethe-Salpeter equation, written schematically as 

T I = K I + tj-0^ JCj(e) A wK (e) Tj(i) , (28) 

where A^k is the two-meson propagator. The diagramatic representation of this equation, 
together with its perturbative solution, are shown in fig. 6. In the case of low-energy hi- 
ll 



















T 




K 


+ 


T 




K 





K 



K 



K 



K 



K 



K 



FIG. 6: (Color online) Bethe-Salpeter equation for the elastic irK amplitude: full equation (top) 
and perturbative solution (bottom). 

teractions, the treatment of the Bethe-Salpeter equation can be enormously simplified, as 
pointed out by Oiler and Oset[14]. The fact that the kernels /Cj involve only effective contact 
interactions and s-channel resonances makes the two-meson propagator A^x to depend just 
on s and eq.(28) can be rewritten as 



T I = K I -}C I Q Tj 

a - ■ d " 



1 



(29) 

(2tt) 4 {( q+ k)/2+£]2 - Ml [{ q+ k)/2-^ ^Mj ' (30) 
The function Q diverges and has to be regularized. As discussed in appendix B, this 
brings into the problem one free parameter for each isospin channel. The corresponding 
regular functions are denoted by Cli and the solutions of eq.(29) become 



(31) 



1 + fi 7 /Cj 

Above threshold, the functions Qi are complex and written as fli = Ri + i I . They involve 
a loop phase loj = tan -1 [//.R/] and their explicit analytic forms are given in appendix B. 

The amplitudes Tj do respect unitarity and, below the first inelastic threshold, can always 
be written as 2 

16/T 



T, 



P 

tan 5i 



sin 5/ e iSl , 

IK,! 



(32) 



(33) 



1 + /Cj ' 

where the real phase shifts 5i incorporate the dynamical content of the interaction. These 
results allow the kernels to be expressed as 

167r tan 5 1 



p (1 + tan<5// tanuj) 



(34) 



The amplitude T is relativistic and we employ the conventions of Refs. [15] and [13]. 
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The unitarization procedure employed in the derivation of Tj generalizes that based on the 
on-shell iteration of the if-matrix [15, 16], which amounts to neglecting the real part of Qj 
and to assuming /C/ ~ [167r/p] tan 5j. 

The behaviors of the amplitudes T and T 3/2 are very different, owing to the presence 
of a resonance in the former. In this case, we use the kernel (27) into eq.(31) and finds 

16tt m K T K 
1/2 P Ml-8-im K T K > [M) 

where the running mass and width are defined by 

Ml = ml + 1 2 R 1/2 , (36) 

m K r K = 7 2 p/(16vr), (37) 

and the free parameter in R was chosen so that M 2 K (m 2 K ) = m 2 .. This yields a unitary 
amplitude T ' which becomes purely imaginary at s = m 2 . Therefore, we call m K the 
nominal kappa mass. Predictions for phase shifts can also be expressed as 

tan^ = ^ (38, 

and the relationship between nominal and running masses is determined by 

o 

Tfl — S 

- 1 + tan<5 1/2 / tancJ 1/2 . (39) 



M 2 -s 1/2/ 1/2 

if-matrix unitarization corresponds to making A4 2 — > m 2 and deviations between both 
approaches are quantified by the factor (1 + tan 5 1/2 / tanc^ 1/2 ). 

Real and imaginary parts of the amplitude T ' eq.(35), together with the corresponding 
ii'-matrix approximation, are given in fig. 7, for the choice m K = 1.2 GeV. In both cases, 
amplitudes become purely imaginary at s = m 2 and peaks occur at lower energies. As one 
discusses in section IX, these shifts in the peaks are a direct consequence of chiral symmetry. 
The figure shows that the ^-matrix approximation is a crude one and that the role played 
by the running mass is important. 

In figs. 8 and 9 we display the isospin dependence of Si and \Tj\ 2 . It is worth noting that, 
by construction, 5 passes through 90° at s = m 2 and, again, the i^-matrix approximation 
is crude. In the case of the 1 = 3/2 channel, the free parameter in R. i/2 was fixed by 
imposing that the predicted phase shifts around s = 0.9 GeV agree roughly with those given 
in Ref.[17]. This gives rise to huge differences between full and ^-matrix results. 
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FIG. 7: (Color online) Real (full lines) and imaginary (dashed lines) components of the amplitude 
T obtained by means of eq.(35) (blue) and in the .fT-matrix approximation (red). 




s(GeV 2 ) s(GeV 2 ) 

FIG. 8: (Color online) Predicted phase shifts for T 1/2 and T 3/2 amplitudes, together with iT-matrix 
results. 

In the decay D + — ► K~7T + 7r + , elastic amplitudes Tj contribute to final state interactions 
only. They are always accompanied by the two-body propagator Cli, as indicated in fig. 2 
[please see also eqs. (44-47)]. Real and imaginary parts of the products CljTj are displayed 
in fig. 10. In order to clarify the meaning of this figure, one notes that, in the iT-matrix 
approximation, Cl — > —ip/(16n), and [5R H T, 3 H T] is given by [p/(167r) (3? T, — 3? T)]. 
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FIG. 9: (Color online) Dependence of |T 1/2 | 2 and \T 3 | 2 on s, together with the K-matrix approx- 
imation. 

When the real part of Q is turned on, a shift in the curves occur and the phase of flj Tj 
becomes Si+ui. The magnitudes of the Uj may be inferred by noting that, in the i^-matix 
case, 9ftf2T = Q ; T = Oat threshold. Loop phases are, therefore, explicit ingredients of 
the decay amplitude. With future purposes in mind, one notes that the condition —1 < 
f2 T 1/2 < 1 holds for both the real and imaginary components of this quantity. 
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FIG. 10: (Color online) Real (continuous line) and imaginary (dashed line) components of the 
functions 0/ Tj. 
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VI. FSI: PRODUCTION AMPLITUDE 



K 



K" 



FIG. 11: (Color online) Resonance (full line) propagation and decay into ir + K (dashed lines); T 
is the unitary 1 = 1/2 scattering amplitude. 

We consider the possibility, shown in fig.2f, that the resonance can be directly produced 
at the weak vertex. The interpolation between the decay and the observed ir + K~ final state 
is described by the subset of diagrams shown in fig. 11, denoted by U 1/2 and referred to as 
production subamplitude. It involves both the bare k propagator and the unitarized elastic 
scattering amplitude T , given by eq.(35). Reading the diagrams, one finds 
VS/F 2 



ill 



1/2 



s — mi 



[c d (s-M 2 -M 2 K ) + c m (4M|+5M^)/6] 



"l/2 J l/2 



(40) 



This function is complex, owing to the factor [1 — fi 1/2 T 1/2 ], and can be cast in various fully 
equivalent forms[18], namely 



in 



cos 5 



1/2 



1/2 



mf. — s 
cos 5 



1 + 



tan 5 



1/2 



tanu; 



e 1 / 2 



V2 g« 1/a 



1/2 

sin 5 



2 e i5 ^ , 



(41) 
(42) 



Ml-s 

g = -(V3/F 2 ) [c d {s-Ml-M 2 K ) + c m (4M 2 K + 5M 2 )/Q] 

Results have very little model dependence and show that, as predicted by Watson's 
theorem, the phase in the production amplitude is 8 1/2 , the same as in the elastic process. 
On the other hand, the magnitude of Il 1/2 is determined by both 5 and u 1/2 . It is important 
to stress that the full equivalence among these different forms holds only for the running 
mass and width given by eqs.(36) and (37). If other forms for these functions are employed, 
consistency is lost and results are no longer under control. As in the elastic case, the function 
n i/2 contains the A^-matrix approximation as a particular case. In Fig. 12 we show the 
behavior of the function iU 1/2 (s), together with the corresponding A- matrix approximation, 
and notes that the differences between both sets of results are important. The figure also 
indicates that bumps in the former are more pronounced. 
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FIG. 12: (Color online) Modulus (full line), real and imaginary components (dashed and dot-dashed 
lines) of the function iU 1/2 (s), eqs. (40-42), together with the K-matrix approximation. 



VII. AMPLITUDE FOR D~ 



TT + 7T + K 



We display here individual contributions to the D + decay width from the diagrams of 
fig. 2 , using both the scattering and production amplitudes derived previously. They are 
covariant and expressed in terms of the invariant masses fiij, defined in appendix A. The sum 
A = [A a H — • + Af] is designed to replace the background plus S- wave factor [oq e<*° A + 
af e l ^i Af] mentioned in the introduction, 
diagram (a): 



A a {fi 



2 ; 

7T7V ) 



1 



3V2 



[S a ] Gp cos 2 6 C /i 



2 . 

7T7T ' 



diagram (b): 



A b (fi 2 



K 71 



3V2 



[S a ] G F cos 2 6 C 



{P-q+Ml 



{P-g-Mp [Ml-MD 
M 2 D +M%-2P-q 



x 



2 2 1 

3 ^1/2 ^1/2 (/%7r') + 3 ^3/2 ^3/2(^71 



(43) 



(44) 
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diagrams (c+d): 



Ac+d(fJ>Kir) 




[5 a ]G F cos 2 6 C {P-q'+Ml) 



(P-q'-M 2 )(M 2 K -M 2 y 
M 2 D +Ml-2P-q' 



^1/2 ^i/ 2 (/ i X7r) ; 



(45) 



diagram (e): 



Ml*kn) = ~\ N G F cos 2 9 c (M 2 D -3P-q>) 



(Ml-P-q') (j^-jg) ' 
M 2 D +M%-2P-q' 



X ^1/2 ^i/ 2 (/ i i"7r) ^3/2 ^3/2 (/^fOr) 



(46) 



diagram (f): 



m _ _ — m ^ 



x [ Q ( /U L-M 2 -M|)+c m (4M| + 5M 2 )/6] 1 - T^O&J . (47) 



The process L> + — > X~7r + 7r + is Cabibbo allowed and, as expected, amplitudes share 
the factor G F cos 2 9 c . As far as phases are concerned, one finds three kinds of structures. 
The amplitude A a comes from a tree diagram and is necessarily real. The phases of the 
amplitudes Ab, A c+ d and A e , on the other hand, are contained in the products O/Tj and 
given by (5i + ujj). Finally, as discussed in [18], the phase of Af is 5 1/2 , the same of free 
scattering. 

In the evaluation of Dalitz plots, the amplitude A has to be symmetrized with respect 
to the variables fi 2 Kn and [i 2 Kn ,, since outgoing pions are identical. Allowed values for these 
invariant masses lie in the interval 0.40 GeV 2 < ^ 2 KlT) l^ n > < 2.99 GeV 2 , whereas the diagonal 
of Dalitz plot corresponds to 0.94 GeV 2 < ji 2 KlT = n 2 KlT , < 1.85 GeV 2 . As discussed in section 
V, we assume our results to be valid for s < 1.6 GeV 2 in the two-body channel. 

Isospin 3/2 contributions are present in Ab and A e only. Their dependence on isospin is 
displayed in figs. 13, where it is possible to see that, for I = 3/2, just the real part is relevant. 

Real and imaginary components of the amplitudes At, as functions of ^ 2 KlTl are shown in 
fig. 14, for the choice m K = 1.2 GeV and 5 a = 5b = S c = 1. As the amplitude A a depends on 
fi 2 n , we rewrote it as 



A a (x) = [5 a ] G F co S 2 9 c (M 2 D +2M 2 +M 2 K -x) 



(48) 
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FIG. 13: (Color online) Full (black), / = 1/2 (blue) and I = 3/2 (red) contributions to the real 
(continuous line) and imaginary (dashed line) components of the amplitudes Ab and A e ] the vertical 
scale has to be multiplied by the weak factor Gf cos 2 9c- 



and just A a (^ 2 K7T ) was included in the figure. All contributions have comparable magnitudes 
in this range, with a dominance of diagrams (a) and (f) in fig. 2, which represent the non- 
resonant background and the direct production of the resonance at the weak vertex. We 
note, however, that the latter is rather sensitive to the resonance coupling constants Cd and 
c m and recall that the values adopted here are just illustrative. 



VIII. DALITZ PLOTS 



In order to produce a feeling for the Ai given by eqs. (43-47), we display here their predic- 
tions for Dalitz plots. The plotted quantity is \Ai(jiKir) + Ai(nKn') \ 2 and individual contri- 
butions correspond to diagrams in fig. 2. Our description for the amplitudes Ai is valid for 
invariant masses below 1.6 GeV 2 . As, in the plots, this condition must hold simultaneously 
for both amplitudes, a reliable region around the lower side of the diagonal is selected. 

Diagram (2. a) gives rise to fig. 15. It describes the non-resonating background and one 
learns that it is not evenly distributed along the plot, as sometimes assumed in the literature. 
We show, in fig. 16, individual contributions from diagrams (2.b-f), which involve final state 
interactions. Amplitudes Ab and A c +d share same weak vertices, but the latter is based on 
just the isospin 1/2 kernel, whereas the former contains an admixture of isospins. However, 
their Dalitz plots are very similar, indicating that isospin 3/2 contributions are small. The 
plot corresponding to A e is very different from the other ones because its weak vertices 
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FIG. 14: (Color online) Full (black) and partial contributions to the real (continuous line) and 
imaginary (dashed line) components of the amplitudes Ai\ the vertical scale has to be multiplied 
by the weak factor Gf cos 2 6c- 



are of the vector type and the W is coupled to two mesons. The direct production of the 
resonance at the weak vertex is associated with Af and produces a plot with a rather broad 
peak around 1.3 GeV 2 . Finally, the sum [A a + • ■ ■ + Af] is given in fig. 17, where a typical 
interference pattern can be noted. 



IX. SUMMARY AND CONCLUSIONS 

The low-energy components of the amplitude K~n + 'K + are studied in the frame- 

work of a rather conservative SU(3) x SU(3) chiral effective theory and special attention is 
paid to the resonance k. For practical reasons, the derivation of weak vertices is performed 
using the group S77(4), but without any commitment with the corresponding symmetry. 
In dealing with final state interactions, proper three-meson processes are neglected and we 
remain within the quasi-two body approximation. Our main results are summarized by 
eqs. (43-47), which represent individual contributions from the diagrams in fig. 2. At low- 
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FIG. 15: Contribution from diagram (a), which represents the non-resonant background. 

energies, the decay amplitude is represented by A = [A a + • • • + Af\ and symmetrization 
with respect to final pions is required. Conclusions are presented in the sequence. 

1. degrees of freedom: The amplitudes Ai contain both fixed and adjustable pa- 
rameters. The former class encompasses pseudoscalar masses Mj, their decay constant 
F = F n = 0.093 GeV and the weak constants G F = 1.166 x 10~ 5 GeV~ 2 and cos6 c = 0.9745 
[12]. Adjustable parameters involve resonance masses and coupling constants. In principle, 
three scalar states should be considered, but two of them appear just in the t-channel back- 
ground. As discussed in appendix C, contributions from t- and w-channel interactions are 
very small and can be safely neglected. In this approximation, the only free parameters in 
our results are m K , the nnK coupling constants Cd and c m , and three SU (4) breaking factors 

2. background: The non-resonating background is represented by A a , which is a real func- 
tion. It important to note that this does not happen by chance. The process shown in fig. 2a 
is the simplest possible in D + decays and given by a tree-level diagram. In field theory, tree 
diagrams are real and imaginary components are produced by loops. Of course, it would 
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FIG. 16: Contributions from processes involving FSIs from the following diagrams: (b) top-left; 
(c+d) top-right; (e) bottom-left; (f) bottom-right. 



be possible to dress the primary weak vertex with mesonic loops, but this would amount 
to including higher order chiral corrections. Therefore, at low-energies, the background is 
necessarily real and should not be represented by trial functions of the form [a Aq\. 
Inspecting figs. 14 and 15, one learns that the magnitude of A a is comparable to other 
contributions and that its distribution over the Dalitz plot is not uniform, as sometimes 
assumed. 

3. phases: Our calculation begins with a lagrangian, which yields real vertices only. Loops 
are introduced in a subsequent step, and only then amplitudes become complex. This 
construction process is systematic and one has full control over all imaginary terms and 
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FIG. 17: Contributions from diagrams (a H f). 



understand clearly their dynamical origins. Complex amplitudes are due to final state in- 
teractions and encoded into the functions 0/ and Tj, representing respectively two-meson 
propagators and elastic Kir scattering. As one deals with two isospin channels, in principle, 
four phases need to be considered. However, the 1 = 3/2 channel is repulsive and the corre- 
sponding phases are small. One is then left with just uj 1/2 and 5 1/2 . Both of them are present 
in the phases of Ab, A c+ d and A e , which are identical and given by (5 + oj 1/2 ). Diagram 
2f represents the direct production of the K-resonance at the weak vertex and involves the 
subamplitude II 1/2 , represented in fig. 11. As pointed out in ref. [18], this structure gives rise 
to the same phase as in T 1/2 , namely 5 1/2 . Finally, one notes that another phase is produced 
when the complex I = 1/2 amplitudes are added to the almost real 1 = 3/2 conterparts. 
In summary, the low-energy decay amplitude contains several different energy- dependent 
phases and cannot be well represented by a trial function such as [af e 1 ^ Af] . 

4. two-body isospin channels: Our numerical results indicate that final state interac- 
tions are very important in D + — > K~it + it + . They are present in five of the diagrams shown 
in fig. 2. Three of them involve just the isospin 1/2 channel, whereas the other two depend 
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on both components. These cases were studied in fig. 13 and one notes that isospin 3/2 is 
relatively important in A e only. 

5. i^-matrix approximation: This problem is addressed in section V. The main advan- 
tage of the i^-matrix is its simplicity. Our expressions encompass the i^-matrix approxi- 
mation, since it amounts to neglecting the real part of the two-loop propagator and setting 
uji = 7r/2. In the isospin 3/2 channel, this procedure is disastrous. For / = 1/2, on the other 
hand, it gives rise to reasonable qualitative predictions. One should bear in mind, however, 
that the phase u 1/2 influences final result in two different ways, since it both helps shaping 
the K7r amplitude and enters directly into the expressions for Ab, • • • , Af. 

6. dynamics: We assume our results to be reliable for Kn invariant masses between 
0.4 GeV 2 and 1.6 GeV 2 . The magnitudes of the Ai are comparable in this range, ss shown 
in fig. 14. There, it is also possible to see that the non- resonant background and the direct 
resonance production dominate. However, numerical results are particularly sensitive to the 
coupling constant q and, for the time being, we take this conclusion as provisional. 

7. Breit-Wigner expressions: Our results involve functions which are akin to the usual 
Breit-Wigner ones. However, as we discuss in the sequence, differences between them are 
very important. These functions are hidden in the two-body amplitudes T 1/2 , present in 
diagrams (2 b, ... ,e) , and also contribute to the Tl 1/2 in (f). Eq.(35) can be rewritten as 

7 2 

Tl/2 = i(ml-s) +1 iR 1/2 ]- i[ 7 2 p/(167T)] ' (49) 

where R 1/2 represents off-shell effects in the two-meson propagator. The function 7 2 is given 
by eq.(27) and can be rewritten as 

1 2 = 3h 2 /F 2 + [aA + c d p}(m 2 K -s) , (50) 
h=[c d (m 2 K -M 2 -M 2 K )+c m (4M| + 5M 2 )/6] , 
A = (1/4F 2 ) [(4-3p 2 /2) S -4 (M 2 + M^)] , 

(3=[3c d (m 2 K -s)-6h]/F 3 . (51) 

In this expression, h 2 is a 0(q A ) effective coupling constant, A is the 0(q 2 ) leading term 
in the chiral amplitude and (3 is a 0(q 4 ) background. A parameter a has been introduced, 
so that the leading chiral contribution could be turned on or off. The usual Breit-Wigner 
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expression can be recoverd from eq.(49), by going to the .fT-matrix approximation (R 1/2 — > 0) 
and by choosing (a = 0,Cd = 0). This yields the curve BW in fig. 18, with its well known 
shape. The reintroduction of R 1/2 produces no visible effects. The choices (a = 1, q = 0) 
and (a = 0, q = 0.032 GeV) give rise respectively to curves a and (3. In both cases, one 
notes huge enhancements in the region 0.8 GeV 2 < s < m 2 . At threshold, on the other 
hand, the chiral hierarchy is respected, since a — > 0(q 2 ) and (3 — > 0(q 4 ). Finally, the curve 
xBW is produced by (a = l,Cd = 0.032 GeV) and is directly related with those discussed 
in section V. 




s(GeV') 



FIG. 18: (Color online) Relationship between a usual Breit-Wigner (BW, dashed line) and its chiral 
generalization (\BW , continuous line); the other curves correspond to the choices a : (a = l,Cd = 
0) and p : (a = 0,c d = 0.032 GeV) in eq.(50). 



8. chiral symmetry: The phase shifts predicted by eq.(49) pass through tt/2 at s = m 2 .. 
This feature defines the nominal mass of the resonance and is completely independent of 
parameters adopted. On the other hand, the implementation of chiral theorems requires that 
amplitudes be represented by polynomials which have well known values at threshold and 
grow close by. The form of the curve xBW in fig- 18 corresponds to a compromise between 
those two features. It is important to note that this kind of behavior cannot be obtained 
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by adding polynomials to usual Breit-Wigner expressions. Rather, it derives directly from 
the unitarization of contact interactions added to resonance poles, as discussed in sections 
IV and V. A study of pole movements induced by this procedure is in progress and will be 
presented elsewhere. 
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APPENDIX A: KINEMATICS 

Initial and final momenta are uniformly represented by capital and low-case letters, 
two-body system: In the description of the two-body reaction ir(Q) K(K) — > 7r(g) K(k), 
Mandelstam variables are defined by s = (Q + K) 2 = (q + k) 2 , t = (Q — q) 2 = (K — k) 2 , 
u = (Q — k) 2 = (K — q) 2 and satisfy the condition s+t + u = 2 (M]^ + M 2 ). In the center of 
mass, one has 



t = ~{sp 2 /2) (l-cos0) , 



(Al) 



u = (M 2 K -M 2 ) 2 /s - (sp 2 /2) (l + cos0) , 



(A2) 




(A3) 



In performing S'-wave projections, one uses 




(A4) 



and finds 



t^t = -sp 2 /2 , 



u^u = (M 2 K -M 2 ) 2 /s-(sp 2 /2) , 



(A5) 



[l/{x-m 2 



)] 



o 



(1/sp 2 ) \n[l + sp 2 /(m 2 -x -sp 2 /2)], 



(A6) 
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with x = t,u. 

three-body system: In the D + decay, momentum variables are defined as D + (P) — > 
7T + (g) 7r + (g') K~{k). The W propagator splits diagrams into two parts and the momentum 
(q) is associated with the pion in the same sector as the D + . Invariant masses are described 
by 

a4 = {q+q')\ A, = (v+k) 2 , = (q'+k) 2 , (A7) 

and one has fj, 2 n + [i 2 Kn + fi 2 Kn , = M|, + 2 M 2 + M\. Scalar products are then given by 
P-q=(M 2 D + M 2 - ^,)/2 and P ■ q> = (M 2 D + M 2 - ^)/2 . 



APPENDIX B: TWO-MESON PROPAGATOR 
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FIG. 19: (Color online) Two-meson propagator. 



The two-meson propagator is shown in fig. 19, for a system with total momentum X 
Q + K = q+k. We also use the combination £ = (Q—K)/2 and define 

[1; £"] 



[Ikk; I% k \ - J 



(2tt) 2 [{£+X/2) 2 -M 2 } [{l-X/2) 2 -M\ 



Bl 



In terms of Feynman parameters, these integrals read 



[I-kK] I% k ] 



(4tt) 



n 00 ■ — (n w - rr 01 

LL irKl 2 v * K irK) 



da \-(l-a) m } [-a n ] ]n[D vK /A 2 ] + 



D nK = (l-a)Ml + aM 2 K -a{l-a) X 2 



(B2) 



where the ellipsis indicate an infinite quantity associated with dimensional regularization. 
Multiplying I^ K by X^ in eq.(Bl) and manipulating the integrand, one finds the following 
useful result 

. . \r 2 - \ri 

(B3) 



M 2 — \ I - 



j£2 "7r_fsT 

The two-meson propagator given by eq.(30) is written as Q = — (1/167T 2 ) [L(s)+Aoo], 
where the function L(s) is explicitly given below and Aoo is an infinite constant that has to 
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be removed by renormalization. In this procedure, the function Q is replaced by 

flj = -(l/167r 2 )[L(s)+c,] , (B4) 

where the cj are constants which depend on the isospin channel. They are chosen by tuning 
the predicted phase shifts Si(s), eq.(33), to experimental results at a given point s = sj and 
one imposes Si(sj) = 5j Xp (s/). When a resonance is present, a rather convenient choice[18] 
for si is the point at which the experimental phase is n/2. 
The function L(s) entering eq.(B4) is given by 



for s < (M K - M n ) 2 : 
~a- 1 



Lis) 



-p log 



7/ 



for (M K - M n ) 2 < s < (M K + M n ) 2 : 

L(s) = p [tan -1 a — 7r/2] — r] , 

for s > (M K + M n ) 2 : 
"1 -a 



L(s) = p(s) log 



l + <7 



- 1] + t7ip , 



o- = ^\s-{M K + M n y\/\s-{M K -M n y\ , 

77 = 2 - p4 - M w 2 )/ S ] loglM^/M,)] . 

These results allow the renormalized two-loop propagator to be written as 

1 



(B5) 
(B6) 

(B7) 

(B8) 
(B9) 



Q 7 = 



16tt 2 



{^t[L{s) - Cl }+i^L{s)} 



= R^+iOls-iMK + M^) 2 } I(s) 



(BIO) 



The imaginary component is very simple and reads I(s) = — p/(16n). We define a loop 
phase uji by 

tanu;/ = I/Rj . (Bll) 

For the channel / = 1/2, the constant Ci/ 2 is chosen so that 

1 



16tt 2 



■ [L(s)-L(mD] 



(B12) 



and, by construction, R 1/2 {m\) = 0. In the / = 3/2 channel, fit to data requires c 3/2 
140 GeV. The functions Ri, I and ui are shown in fig. 20. 
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FIG. 20: (Color online) Left: real (continuous line) and imaginary (dashed line) components of the 
two-meson propagator; right: propagator phases. 

APPENDIX C: BACKGROUND AMPLITUDES 

The resonance-exchange amplitudes calculated in section IV have the general form 



It = ~ ^4 l^ m l + C ") + Ck)} + c 2 d (x+ml) + c d (C»+C£)} , (CI) 
where x — t,u, and 

C: = Cl = —c d (M 2 + M 2 K ) + c m (4M 2 K + 5M 2 )/6 , (C2) 

C° = -2 (c d -c m ) Ml, C° K = -2 (c d -c m ) M 2 K , (C3) 

= -2c,M2 - c m (2M 2 K -llM 2 )/6 , (C4) 

C| = -2c d M 2 K + c m (WM 2 K -M 2 )/Q . (C5) 

Using the results for 5-wave projection given in appendix A, the background amplitudes 
entering the kernels of section IV are written as 

B l/2 = 3 S;/4 + B° - Bf/6 - B-JA , (C6) 

B 3/2 = B» - B*/Q + B*/2 , (C7) 

B: = - ± [c 2 d (s+m 2 K ) + c d (C:+C" K )} , (C8) 



+c 2 d (x +m 2 a ) + c d (C: + C a K )} . (C9) 
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FIG. 21: Background amplitudes. 
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